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ABSTRACT 


FLUID CONTROL SYSTEM BEHAVIOR IN A RANDOMLY VIBRATING ENVIRONMENT 
LCDR. Franklin F. Alvarez, U.S.N. 


Submitted to the Department of Mechanical Engineering on August 15, 
1966 in partial fulfiliment of the requirements for the degree of 


Doctor of Science. 


A fluid control system is analysed while operating in a 
vibrating environment. The introduction and transmission of noise 
signals in the way of unwanted pressures are examined for passive 
transmission line elements and an active flapper-nozzle control valve 
when subjected to environmental random or periodic translational 
vibrations. Methods are presented for predicting the levels and 
frequencies of induced pressures using distributed and lumped parameter 
techniques, including some damping effects. Analytic results evalu- 
ated with the help of a digital and analog computer are presented 
along with experimental data obtained for vibrating pneumatic trans— 
mission lines and for a small pneumatic control valve. Also included 
is a procedure for testing a given system's susceptibility to vibra- 
tion induced noise sources. 

It is found that pressure peaks proportional to the fluid 
density, speed of sound in the fluid, and degree of excitation are 
generated in transmission lines over a wide band of frequencies. The 
peaks occur at the system resonant frequencies, and are due to induced 
fluid flow caused by pipe boundary nonuniformities (bends and area 
changes) moving relative to the fluid. The noise signals are excited 
by either longitudinal or transverse piping vibrations making the 
system susceptible to noise introduction regardless of its orienta- 
tion with respect to the vibrating field. 

Noise is introduced into the control valve by means of the 
vibrating supply line, inertial forcing of the valve flapper, and for 
high pressure systems, by vibrations of the transmission line to the 
load. The greatest noise transmission takes place at low frequencies 
and near the flapper resonance. Force feedback by means of a bellows 
arrangement decreases the amplitudes at the lower frequencies, but 
reinforces those near the flapper resonance and can cause noise power 
amplification. 


Thesis Supervisor: Henry M. Paynter 


Title: Professor of Mechanical Engineering 
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1. INTRODUCTION 


— 


BACKGROUND 


| 


High performance hydraulic and pneumatic control systems are in 
wide use today because of the high speeds of response they offer. 
These systems can be found in missiles, aircraft, and ships, each of 
which is at times subject to random vibrations from either internal or 
external sources such as rocket combustion, engine vibration, or ran- 
dom seas. A fluid control system operating in this environment is 
subject to noise, unwanted pressure signals, excited by induced motion 
of the system structure and of the working fluid itself. How this 
noise is introduced into the system and how it effects the transmission 
of command signals are the two main questions this thesis will try to 
answer. Other important considerations are the dependence of noise 
generation on the orientation of the system to the vibration field, 
coupling between the fluid and structure, and transmission coupling 
into directions other than that of the primary vibrations. 

Studies of fluid control valve forces both steady-state and 
transient have been conducted in the past Ire 2, 3] . However, none 
of the cases studied have included any environmental vibrations of the 
valve assembly itself with inertial forces being of significant impor- 
tance. 

Likewise, transmission line behavior has been investigated [4, Би 6 | 
for a line that is in an inertial system. D'Souza and Oldenburger [7 | 
included some coupling between the piping and fluid, but they still did 


not consider any vibrating environment induced motion. 
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Silberman [8 | did vibrate a tube with fluid flowing through it, 
but this investigation was solely to find the effect on pressure loss 
between two points, and it did not look at any pressures induced by 
the vibrations. 

In Lipner and Fay [9 | one finds the first mention of environmental 
induced pressures. This paper discusses the transient behavior of 
acoustic waves generated by a seismic shock. The simplified analysis 
begins with a pipe of infinite extent and ignores end conditions even 
when boundaries are established. Thus fluid resonances are not 
considered. 

Random vibration studies have been conducted investigating parts 
failures in electronics equipment [10], Structural response to random 
forcing Le, 12 13 | , and fatigue failure of materials [12, 13, 14 | : 
Once again, hovever, there is a lack of information dealing vith fluid 


power elements while operating in a randomly vibrating environment. 


1.2 OBJECTIVES AND PROCEDURE 


The objectives of this thesis are three fold. The first is to 
investigate the sources of noise production and its transmission in 
fluid systems subject to environmental vibrations. The second is to 
analyse the behavior of an active element, a control valve, when it is 
part of this system, and to be able to give procedures for estimating 
its susceptibility to this type of noise. The third is to be able to 
recommend the most profitable avenues for further investigation in this 


area. 
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In order to conduct meaningful experiments in conjunction with the 
analysis, the following procedures were adopted. Air was used for the 
working fluid because of the large quantities of pressurized air 
available and the simplicity of exhausting directly to the atmosphere. 

A nozzle-flapper control valve was decided upon as the active element 

due to its ease of construction and the availability of previous 
analytical information. Also, it was decided to use a choked orifice, 

a pure resistance, for the load as no further insight would be gained 

by having the transmission line terminate in a more complicated impedence. 
Both analog and digital computers were used to help interpret the 
analytic results. 

Throughout this thesis the units of pounds mass, pounds force, 
inches, and seconds are used for mass, force, length, and time respective- 
ly. This introduces in Newton's Second Law the factor go? relating pounds 
force to pounds mass, which is equal to 386 lbm in per lbf pecan 

The thesis contains in Chapter 2 a study of vibration induced 
noises in transmission lines subject to both longitudinal and transverse 
vibrations. The last section of the chapter contains a summary due to 
the relative importance of this chapter to the thesis. Chapter 3 
includes the analysis of a control valve operating in a vibrating 
environment, discussing the sources and transmission of noise. Also 
included is a procedure for testing a given system's susceptibility to 
vibration induced noises. The final conclusions and recommendations are 
presented in Chapters 4 and 5 respectively, and supplementary work is 


contained in the Appendix, a 
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2. VIBRATION INDUCED TRANSMISSION LINE NOISE 
2.1 INTRODUCTION 


The basic system to be investigated consists of a supply line 
leading from a constant pressure reservoir to a flapper-nozzle 
control valve. The load controlled by the valve, in this case a pure 
resistance, is located at the end of a transmission line, Fig. 2.1. 
The entire system is assumed to be on a rigid base that can experience 
translational vibrations in all directions. Since the base is rigid, 
all supports, where joined to the base, experience the same motion 
independent of positions. The entire system is subject to induced 
vibrations due to the forced motion of the base that for this inves- 
tigation assumes the role of a rocket frame, plane fuselage, or ship 
hull excited by periodic or random accelerations. Due to the induced 
vibrations of the structure, noise is introduced into the system by 
means of motion of the control valve flapper relative to the nozzle 
and by the motion of the structure, containing the fluid, relative 
to the fluid. The valve flapper and actuator comprise a second order 
system and therefore are susceptible to noise when excited near the 
resonant frequency. Since their behavior is dependent upon the control 
valve properties, this source of induced noise will be investigated 
in Chapter 3. The rest of this Chapter will be devoted to fluid 


pressures induced by structure motion. 


222 CONCEPT 
When a vibrating pipe contains a fluid under pressure, the 


boundary condition at the interfaces is that the fluid has the same 
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motion as the piping. Thus for pipes experiencing vibrations along 
their longitudinal axis, if the contained fluid is of low viscosity 
such as air, the boundary layer is thin, and the main mass of fluid 
is at rest. However, the fluid at end junctions such as bends or 
area changes would be experiencing motion because of the longitudinal 
vibration of the solid boundaries. These junctions act as flow 
sources to the fluid in the same way as a piston oscillating ina 
stationary pipe. These flow sources excite fluid vibrations causing 
quite high pressures to be obtained at the resonant frequencies. 

For liquids, the boundary layer is thicker than for air, but the 
Same general behavior takes place. Cavitation bubbles can form at 
the solid boundaries, however, reducing the effectiveness of the 
structure motion. 

When the induced vibrations are transverse, all the fluid is 
forced in a transverse motion along with the pipe. However, due to 
mechanical coupling the pipe sets up longitudinal vibrations at the 
bending resonant frequencies and at double the resonant frequencies. 
Then as before, any junctions act as flow sources to the fluid and 
impose longitudinal vibrations on top of the transverse motion. 

For this investigation, the base is assumed rigid and all support 
excitation is the same independent of position. However, in actual 
practice the piping structure can be excited independently at different 
locations especially if the distance between points is large. Thus 
each junction can act as an independent flow source, and the total 
fluid response will be the superposition of the responses to each of 


the sources. Likewise, if the excitation is both longitudinal and 
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transverse, the pressure response is the superposition of the two 
types, and each can therefore be treated separately in the analysis. 


These concepts will be expanded upon in the following sections. 


2.3 METHOD OF MODELLING VIBRATION INDUCED FLOW SOURCES 


—< s — mm me ee ee  — ee Be A 


The representation of the flow source concept by a physical 


system which can be modelled is shown in the following sketch. 


= 








+ 
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The representation allows a vibrating bend or area change to be 
modelled as two transmission lines forced at their junction by a 

flow source. Thus normal analytic techniques such as developed by 
Paynter [15 | for energy transmission elements can be used. The 
forcing functions are then obtained when the additional conditions 

for the interaction between the fluid and piping are included in the 
analysis. These forcing terms which are peculiar to systems under the 
influence of environmental vibrations will be identified in a sample 


problem formulated in the next section. 
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2.4 DEVELOPMENT FOR LONGITUDINAL VIBRATIONS 
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First to be considered is piping forced in a longitudinal 
motion. Since the fluid is forced at junctions it is the conditions 
at these junctions that are of most interest. Thus the analysis 
leads to the use of transmission matrices between junctions which 
can then be cascaded to obtain the overall effects. This method will 
be used here. A development of uniform, one-dimensional, linear, 
distributed system transmission matrices, along with extensions to 
systems in noninertial frames and with mean transport, is contained in 
Appendix A. 

In order to apply the transmission matrix, Eq. (A.1-3), to the fluid 
in a transmission line it is first necessary to obtain a relation for 


the variables in differential equation form, That is, 


dp _ _ 

DE ZW 
and 

dw _ 

dx Kaes 


This is necessary in order to identify the system parameters Z 
and Y when the power variables p, pressure, and w, mass flow rate, 
are used. 

Friction and damping will be ignored for the present as well as 
viscous transfer of motion from the pipe wall to the fluid. 

Using the continuity equation for a small elemental control volume 


of length dx in a uniform pipe gives, Fig. 2.2, 


0) — 77 7: 


0X öt 
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or 


ow др ӘА (50 
A s Tp ot 


where p is the fluid density and À is the flow area. 
When using the acoustic assumptions of small pressure fluctuations 


and an adiabatic process, the equation of state becomes 
dp = E do 


where ß is the bulk modulus of the fluid. 
The change in area is found by conducting a stress analysis on 


a section of the thin-walled elastic pipe. This gives 


o(D_- D,)dx = D.pdx 
or 


do(D_- р,) = D dp 


D 


where o is the circumferential stress in the pipe and D. and D, are 
the outside and inside diameters of the pipe respectively. 


Now the change in circumference with area is 


where E is the modulus of elasticity and C is the circumference. 


Therefore the change of area with pressure becomes 





AJD: 


2A 
dA = D. dp. 
EG” - 1) 
i 


Putting this and the equation of state into Fq. (2.4-1) gives 


Św _ ФА Әр S 
x 8 ət (2.4-2) 
e 
where 
E 1 2 
EE = (2.4-3) 
Bo B E(D /D; 1) 


It is to be noted that any compliance of the pipe is taken into 
account in Eq. (2.4-3) by decreasing the fluid equivalent bulk 
modulus and thus decreasing the speed of sound in the fluid. This 
effect is important in liquids which have a high bulk modulus. 
However, it can be neglected in air unless the pipe is very compliant 
and has a low modulus of elasticity. 

Now applying the momentum equation to the same elemental control 


volume gives, 


Z ds 
EE: g 9 x g dt 
O O 
or 
Ha ow, w Ov 1 ди ди 
дх SA IX S 8 À x n 8 À dt 020: 


where v is the velocity of the fluid particles. 


From continuity, 


śm, N 9p - _ PA Әр 
ox pA дх KE x ` d ot ” 
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mp. ge mg = = on w ww 
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Putting the equation of state in this expression gives, 


aa | ор, у Әр (2.4-5) 


ox Bo ot B, 0X 


Now putting Eqs. (2.4-2) and (2.4-5) into Eq. (2.4-4) and noting 


B 
from the equation of state that Bo 2a = Bo DE с one can obtain 
2 
( NA „3 1) op. mov op 1 ди 
2 dx 2 3t 7” A öt 
© c o 
Since — << 1, this becomes 
opm l aw (2.4-6) 
öx g A öt ` 
o 
Using the differential operator D = ər Eqs. (2.4-2) and (2.4-6) 


can now be written in the desired form. Thus, 





d — – 
Ce Zw 
and 
dw 
where 
D o 
2-—— and Y = 
g A 2 
O c 


Then 2. = İZ/Y, defined as the characteristic impedence, becomes 


c 


Z = — 


C BÊ 


and y = \vz, the propagation operator, becomes 


o İM 


= 
It is also desirable to express these relations in terms of a 


relative flow variable wee Thus from Eqs. (2.4-7) one obtains 
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dp _ _ z 
тр Zw 20 ADX_ 
and (2.4-8) 
dw, 
dk? 
since w= w + po ADx . 
r O O 


x, is the component of motion of the base in the direction of the 
longitudinal axis of the pipe, Po is the mean fluid density, and 
Wë Y, Los and y have the same values as above. 

The two sets of Equations, (2.4-7) and (2.4-8), can now be put 
directly into there corresponding transmission matrix forms, 


Eqs. (A.1-3) and (A.2-3). Thus for the two cases they become 


p (0) coshyx Z „sinhyx p (x) 
w(o) = sinhyx coshyx w(x) 
c 
and 

p(o) coshyx Z sinhyx p(x) Z sinhyx(p ADx ) 

= 1 + 
56 Zu = coshyx w, (x) (coshyx-1) (p _ADx ) : 

(2.4-10) 


If transport is important, the transmission matrices can be 
obtained from Eqs. (2.4-9) and (2.4-10) by using the multiplicative 
factor as discussed in Appendix A, 

For the piping containing the fluid, the variablesf, force, 
and €, absolute displacement will be used. Once more friction and 
damping will be ignored at present. Using Newton's Second Law for an 


element of pipe, Fig. 2.3, gives 


p A 2 
DPP 


~ — dx = dx 
Bo ot 
Or 
p A 2 
of Ben dire (2.4-11) 
EE: g 2 
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where — and En are the pipe cross sectional area and pipe density 
respectively. Now since f is positive when compressive, Hooke's Law 


becomes 


Í 
E m, (2.4-12) 


vhere E is the modulus of elasticity. 
Revriting Eqs. (2.4-11) and (2.4-12) using differential operators, 


one obtains 


df ` 
dx Be 
and (2.4-13) 
dš _ _ 
dx ` 2 
where 
2 
e À D 1 
7, DOD uva , 
g EA 
O P 


Thus, the characteristic impedence and propagation operator are 


p E 
Z = АР |-2- 
ср P = 
and > 
= D + respectively. 
Yo BE P y 


Once again it is desirable to obtain these relations in terms 


of a relative flow variable d Thus Eqs. (2.4-13) become 
-> = = 26, = zx, 


and (2.4-14) 





= aes 


and Z, Y, Zee and in have the same values as above vhen the 
substitution, Ë = z + Xo? is made. 

Then since Eqs. (2.4-13) and (2.4-14) are in the proper forms of 
Appendix A, they can be directly expressed as transmission matrices. 


Thus for the two cases, 


f (o) coshy x Z sinhy x f (x) 
Yp cp Yp 
E(o) == Ы Х coshy x E(x) (2.4-15) 
SR p p 
D 
and 

f (o coshy x Z sinhy x f (x Z sinhy x x 

) D . Le (x) = ıı 

= d + 

O —- sinhy x coshy x x coshy x-1)x 

EL ) m fe > EL ) ( b ) 5 


(2.4-16) 


Using the relations developed above it is now possible to 
investigate the behavior of fluid contained in a vibrating pipe. 
Consider first the simple problem of a pipe open to the atmosphere at 
one end and closed at the other. Let the open end experience vibra- 


tions in the direction of the pipe axis. 


E 


O E 


The end conditions are, 
£ (o) = 0, plo) = 0, f(L) = -p(L)A, and w (L) = p ADE (L). 


Using Eq. (2.4-16) gives 
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A : 
(1 - coshT E SÉ sinhT P (L) 
2 0 2.00 
üə ein f coshf 


where [ =y L. 
P P 


Then from Eq. (2.4-10) one ohtains, 


2. 
() = | <oshreoshr + 3 (p AD) sinhTsinhr | p(L) 
cp 


x 
+ O SD 2 51пһГсовҺГ (р Арх, ) 


or 


- 2 PA sinh F Dx 
p(L) = - 
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7, 
cosh[coshr RE (p də TEKI 
р 2 o p 


where [ = yL. Now if the excitation is measured in terms of a g of 


: . : Z 
acceleration, which is 386 in per sec , this becomes 


-386Z p A sinhrG 
C O 


1-..vv “ — —— ———— 


D | coshTcoshf + © p A”DsinhTsinhr İ 
De o p 


The frequency response of Li can be obtained by substituting iw 


for D. Thus, 


The last two terms of this equation are the flow sources caused 
by the pipe's vibrating. The next to last term is due to the 
closed end of the pipe moving relative to the open end, while the 
last term is due to the translational motion of the pipe. 
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This equation shows that the pressure levels are proportional 





to fluid density, equivalent speed of sound, and excitation. High 
pressures will be induced at the frequencies at which the denominator 
minimums occur. Relative minimums occur at the resonant frequencies 
of both the fluid and the pipe, and absolute minimums occur at the 
overall system resonances. Damping would decrease the induced 
pressure levels somewhat, but would not change the overall appearance 
of the response. 

Another system might have the pipe mounted to the base by means 
of cantilevered supports. The pipe might then be supplied from a 


reservoir and exhaust through a choked orifice. Such a system, Fig. 


p) 


2.4, was analysed, and the theoretical frequency response for 206 








is plotted in Fig. 2.5. It can be seen that relatively high 
pressures have been induced by the vibrating environment at all the 
odd numbered fluid resonances. The large peak near 800 cps is due 
to resonance of the pipe on the cantilever supports. This also 
causes the fluid resonance near 850 cps to be much higher than the 
others as it is reinforced by the near resonant motion. 

When the cantilevers are stiffened by a factor of two, the 
resonant peak is moved far enough to the right that it is off the 
Figure. Also it is noticed that all but the first fluid resonant 
pressures have been decreased slightly. 


The Figure also contains a plot which includes a fluid mean 
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transport velocity of one tenth the speed of sound for the stiffened 
system. It is seen to have very little effect. This indicates that 
small values of transport are not too important. 

It is to be noted that in the Figure the resonant peaks are 
damped by the orifice resistance and by an assumed damping of 


approximately .04. 


2.5 DEVELOPMENT FOR TRANSVERSE VIBRATIONS 
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The generation of pressures in pipes undergoing transverse 
vibrations is more complicated than in the previously developed 
longitudinal case. The pipe does exert a small component of force 
on the fluid in the direction of the pipe axis due to the bowed shape, 


4 


— 


ek — 
— £] = 


> en A ŞB dx = Pon dx (IX ay 
— Acceleralvon a ee x (sx) 2 
X O ót 


but this does not produce any significant pressures as the slope 
(force component) is too small. Coriolis and radial accelerations 
are also insignificant for this particular configuration. Experi- 
mentally, however, relatively high induced pressures were observed 
along with longitudinal vibrations of the pipe ends even though they 
were braced. Significant end motion was observed at the pipe bending 
resonant frequencies and at double these frequencies indicating 
considerable coupling between transverse and longitudinal motion. 

One means of coupling can be attributed to the pipe maintaining 


approximately the same length while undergoing transverse vibrations. 
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The pipe curvature thus causes longitudinal motion of the ends in 
order to maintain this length. This type of coupling vill be 
investigated nov for a pipe under tvo different types of end conditions, 
pinned-pinned and clamped-clamped. In so doing the maximum transverse 
displacement of the pipe will be determined and used to find the 
change in projected length. The Bernoulli-Euler model for beams will 

* 
be used. 


Consider first a pipe pinned at both ends. Jacobsen and Ayre 
[17| give 


Y = Cy COSHX + c_sinux + c_coshux + c,sinhux (2.5-1) 


2 3 4 


for the space part of the solution of the beam equation, 
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where u = ( ) , y(x,t) is the transverse displacement, 


8 ET 
Y(x) is the space dependent part of the displacement, and the c's 
are the integration constants. I is the pipe cross-sectional moment 


of inertia. 


* 
Reference [16] p. 332 states that if the radius of gyration of the 
cross section of the beam, y = VI/A_, is small compared to the 


wavelength investigated, then the BErnoulli-Euler model is adequate. 
kx 
The normal beam equation without any fluid terms can be used as 


Long [18] has shown that the fluid has very little effect on 
transverse vibrations. This is especially true for air. 
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Nov if the platform holding the end supports undergoes motion in 


a transverse direction, the end conditions will be 


2 
at x = 0, Y = Lë and 77 = 0, 
dx 
dv 
ünü at x = L, Y= y_ and™ = О 
O 2 
dx 


where Yes is the transverse motion of the supports. 
Putting these expressions into Eq. (2.5-1) and into its second deriva- 


tive gives 
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Since this type of forcing causes resonances only at the 
symmetric, odd numbers, modes the maximum displacement will always 


occur at mid span. Thus writing Eq. (2.5-1) at x = = in terms of 


relative displacement, 75” Y ə it becomes 


L 1 1 
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Written in terms of g's of acceleration the frequency response becomes 


y. D) 386 1 1 
mə” o = OL - 1 (2.5-2) 
w 2 cos 2 2 cosh 7 


This then gives the maximum relative transverse displacement of the 
pipe at various frequencies. 


Now following the same procedure for a pipe clamped at both 
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ends, reference [17] gives 


Y = c, (cosux + coshux) + c, (cospx — coshux) 


+ c.(sinux + sinhux) + c,(sinux - sinhux) (2.5-3) 


as a more convenient form for the space part of the beam equation 
solution when complicated end conditions are encountered. 
The end conditions are now 


Y 
Be e 0, yl = ve and = = 0, 


dY 
and at x = L, ram and a = 0. 


Putting these into Eq. (2.5-3) and into its first derivative gives, 


no 
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As before, writing Eq. (2.5-3) at x T in terms of relative 


displacement and g's, gives for the frequency response 
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Reference Di States that the shape of a pinned-pinned beam forced 


at the supports at resonance approaches the normal-mode shape, 


e пих e e © 
sin ə corresponding to the resonant frequency. Since it is 


resonance that is of interest the pinned-pinned pipe displacement is 
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taken to be 


y, = y,, sin = sinwt (2.5-5) 


where n = 1,3,5,... is the number of the mode and ye is the relative 
maximum displacement obtained from Eq. (2.5-2) for a given assumed 
damping. 

Using this equation the change in length of pipe with pinned- 


pinned end conditions can be found as follows: 
y 
as ay 
dy 
Ser — 2 
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Now the assumption of this section ís S = L so that 
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However since L#L , with little loss in accuracy this can be written 
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Putting in Eq. (2.5-5) gives 
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Therefore, 
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Now following the same procedure for the clamped-clamped pipe 


the normal-mode shape DI is 


2557: İ sinhux - sinux + a(coshux - cosux) | sinut 


where 


_ SinhuL - sinuL 


— cosuL - cosuL : 


and AE the maximum relative displacement, is obtained from Eq. 
(2.5-4) ° 
Differentiating, squaring, and integrating the expression in Eq. 


(2.5-6) gives 
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AL m ( kr (1/20 /2sinh2uL + uL)- sinuL cosuL + coshuLsinuL) 
+ 1/2(uL + sinuLcosuL)+ 2a [1/4 cosh2uL + 1/2 + 1/4 cos2uL 


+ 1/2(sinhuLsinuL - coshuLcosuL) - 1/2(coshuLcosuLtsinhuLsinuL)| 


+ a” [1/202 sinh2uL - uL)+(coshĻLsinuL - sinhuLcosųL) 


+ 1/2(uL - sinuLcosuL)] 5 (2.5-8) 


Since Eq. (2.5-8) only holds for the resonant frequencies vhich 


are 
ир), = 4.730, 
UL) a = 10.996, 
UL) 5 = 17.2795 

and KL), = 23:302; 

and since cosh 10.0 = sinh 10.0 = 11,013.2 so that a = -1, the Equation 


can be simplified for all the resonant modes above the first. Cancell- 
ing terms it is found that 


y 
AL = IŻ ( eeh L È — 1- 1/2 cos2uL | ae (2.5-9) 


Thus the change in pipe length is determined for the two types of 
end conditions, and Table 2.1 can be constructed using Eqs. (2.5-7), 
(2.5-8), and (2.5-9). 
Now since su: = 1/2 - 1/2 cos2wt, the experimentally found 
end double resonant frequencies can be attributed to this means of 
coupling, and the end accelerations can now be calculated as soon as 
the pipe length and ratio of maximum displacement to length are known. 
Some typical values of 5 vere obtained using Eqs. (2.5-2) 


and (2.5-4) for a sinusoidal platform acceleration of one g and an 


approximate damping ratio of .04, and are presented in Table 2.2. 
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TABLE 2.1 


Pipe Longitudinal End Movement Due to Bending 


Pinned-Pinned Clamped-Clamped 
Mode AL Mode AL 
y y 
1 2347 EL y ee 1 6.35. ( 2 y“ ee 
3 22.21 с 3 57.71 A 
5 61.68 5 144.95 de 
7 120.90 7 271.69 2 
TABLE 2.2 


Ratio of Mid-Pipe Transverse Displacement to Pipe Length 
(pipe is 1/4" copper tubing filled vith air at 200 psia) 





Pinned-Pinned Clamped-Clamped 
Mode Cps L A Mode Cps L = 
1 аг 1 t. 1700100 
3 25 AS x 10° 3 296 27 EST 
5 606 27 2.3 x 107” 5 733 27 1 ән 
7 240 60 ROKKO © 7 276 60 NET 


Now using the results of Tables 2.1 and 2.2 some typical 
values of the ratio of end double frequency accelerations to base 


accelerations are calculated and given in Table 2.3. 
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TABLE 2.3 
Ratio of Pipe End Longitudinal 


Acceleration to Platform Transverse Acceleration 


(pipe is 1/4"" copper tubing filled with air at 200 psia) 





Pinned-Pinned Clamped-Ciamped 
BB cs Acer TE ida. со, әса) =“ 
base base 
1 68 10 1 152 „02 
3 436 1.18 3 592 .34 
5 1212 07 5 1466 „016 
7 480 . 0033 : 552 .00065 


The values from the Table will be a little high due to the 
restraining action of the end supports. However, they show that the 
end accelerations are of a relatively high value for the fírst and 
third modes, but that they drop off quickly after that. This trend 
and the above values will be compared with experimental results in 
Section 2.8. 

While the above analysis gives valid results for the twice 
resonant frequency end motion, only a qualitative explanation can be 
given for the even stronger coupling at the resonant frequency. 

This effect can be due to nonsymmetry in the system since 
an unsymmetrically supported pipe can be forced into longitudinal 


motion by transverse vibrations as shown. 
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Thus any nonsymmetry could cause the longitudinal vibrations 

to begin, and once started, they very probably would be maintained. 
If the platform is excited with components along both the 

longitudinal and transverse axis, the longitudinal end motion will 


be the additive effect due to both. 


2.6 DAMPING AND ITS EFFECTS 

The previous sections have shown that ail the important pheno- 
menon occur at resonances. Therefore, in order to have any feeling 
for the amplitudes involved it is essential to know the amount of 
damping. So far very little has been said about this for several 
reasons. First, most of the mathematical developments become con- 
siderably more complicated when losses are taken into account. 
secondly, it is difficult to determine what form the losses should 
take, i.e. internal friction, external friction, etc. Lastly, the 
correct values for the parameters are very Hard to determine, and 
sometimes vary by as much as a factor of twenty. An example would 
be support damping of beams which changes greatly depending on the 
tigntness of the connection. Thus the most satisfactory way to gain 
an understanding of the damping involved in a given situation is 


experimentally. 
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* 
The definition for lightly damped systems, 


w 
Aw = 260 = — ; (2.6-1) 


n Q 
where Aw is the bandwidth, & is the damping ratio, and el is the 
resonant frequency, gives a means of experimentally determining damp- 
ing by measuring the bandwidths of resonant peaks. This measurement 
is made at the half power point which is ./07 of the peak amplitude. 

A plot of Eq. (2.6-1), Fig. 2.6, shows that as damping increases 
the peak amplitude decreases and the peak broadens out. Brown [19] 
had also shown theoretically that fluid damping is frequency dependent 
and increases significantly as w goes up. 

Eq. (2.6-1) gives good approximate values of damping from 
experimentally obtained amplitude verses frequency plots. However, 
it is desirable to be able to predict the effects of damping on 
amplitude in analytical discussions as well. This can also be done 
by using Eq. (2.6-1) and experimentally obtained values of damping. 
The zero damped infinite peaks of a plot can be rounded off by making 
its shape conform to the correct shape, in Fig. 2.6, for the given 
damping. This will give a good estimation of the peak amplitudes 
if the values of damping used are close to that actually experienced 
by the system. Good repeatable damping values for both the fluid 
and mechanical systems were obtained experimentally as discussed in 


Section 2.8. 


The fluid and mechanical systems in this investigation are all 
lightly damped. 
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2,7 RANDOM EXCITATION 


The developments in sections 2.4 and 2.5 were conducted using 
sinusoidal excitation to the base. Although rockets, airplanes, 
and ships are often subject to narrow-band periodic disturbances, 

a majority of the time the excitation will be of a wide-band random 
nature. A convenient means of characterizing a wide-band random 
signal is to use a stationary process which has a uniform, band- 


limited spectral density. This is known as white noise. 





White noise is a close approximation to the desired physical process, 
and as long as the frequency band of interest of the system to be 
investigated falls between an and Wo it will give the desired results. 

One can develop the relation for a system response spectral 
density when the excitation spectral density is known. A development 
is contained in Crandall and Mark [14] » pages 70-71, and the 


transformation notation used is 


S(w) = 2- Е(тје 1 dt (2.7-1) 


== OO 


where S(w) is the spectral density and R(t) is the autocorrelation 


function. The results given are 


- H(-w)H(w)S_ (w) 
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or 


2 
S (w) = |H(w) | s (a), (2.7-2) 


where is the response spectral density, S (w) is the excitation 
spectral density, and H(w) is the frequency response of the system. 

The response spectral density is relatively simple to obtain from 
Eq. (2.7-2) since the frequency response for the system can be found 
using the techniques developed in Sections 2.4 and 2.5, and S. (ə) is 
simply 


S, 4 < lel so 


O 2 


0 elsewhere 
for the excitation assumed in this investigation. 

It is to be noted that the spectral density covers both positive 
and negative values of frequency and is in the units of quantity 
squared divided by circular frequency. Therefore if experimental 
work is being done and only positive values of frequency are used 
and one is working with cycles instead of radians, then in order for 


the excitations to have the same value a factor of 4m must be used, or 


W. (f) ArS (ə), 


Then 


il 


WCE) Inc) |? W. (f), (2.7-3) 


where və and 7 are the excitation and response spectral densities in 


terms of cps. 


2.8 EXPERIMENTAL RESULTS 


In order to verify the theory developed in the previous sections, 
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experiments were conducted by actually vibrating pipes. The pipes 
were vibrated longitudinally and transversely with both periodic and 
random excitation while containing air under pressures of 100 and 
200 psia. In some cases the air had a small mean transport and in 
others it had no flow. Some preliminary work was done with water in 
the pipes, but this was discontinued due to the difficulties 
encountered in trying to remove the entrapped air bubbles. The 
piping used was standard one quarter inch copper tubing. The vibra- 
tion generator was built from a high fidelity speaker, and the 
pressures and accelerations were obtained by using piezoelectric 
transducers. Data was taken up to 1000 cps since above this fre- 
quency the system response is almost nil. 

Figs. 2,7 and 2,8 show the piping arrangement for longitudinal 
and transverse vibrations respectively. One end of the piping is 
connected by flexible hose to a constant pressure reservoir while the 
other end terminates at the pressure transducer, Fig. 2.9. Additional 
details of the experimental procedures and instrumentation are 
contained in Appendix C. 

Fig. 2.10 shows good agreement between the experimental and 
theoretical frequency responses for a system of total length 51.6 
inches forced by longitudinal vibrations while containing air under 
200 psia pressure. There is a choked orifice near the transducer 
end of .032 inches diameter which gives a mean velocity of 14 feet 
per second or about one one-hundredth the speed of sound. The Figure 


shows fluid pressure peaks at all the odd resonances due to a cosine 
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function characteristic equation, and a peak near 850 cps due to the 
longitudinal resonance of the pipe. The resonances are not multiples 
of each other because of the frequency dependence of boundary layer 
friction which causes the boundary layer to be thinner at higher 
frequencies. The analysis used in this investigation did not take 
this change in boundary layer thickness into account so that a 
frequency correction to the equivalent speed of sound must be made 
from Fig. 5. It is necessary to use this Figure to obtain the 
correct value of c in order to compute the correct values of the 
resonant frequencies. 

From Fig. 2.10 it is seen that high enough pressures are 
generated per g of acceleration at the first several resonances to 
be of concern if the system were subjected to any substantial 
environmental excitation. Also, from the tirst two peaks, using 
Eq. (2.6-1), the fluid damping ratios are found to be approximately 
.03 and .04 respectively. 

Fig. 2.12 is the frequency response for the pipe longitudinal 
accelerations at the nonforced end. This Figure shows good agreement 
in resonant frequency between the theory and experimental results, 
but the agreement in bandwidth is not too good due to limitations in 
the equipment resolution capabilities. The resolution at high 
frequencies is only one tenth as good as at low frequencies. The pipe 


end had to be mass loaded in order to bring the first resonance down 





This Figure was taken from Brown 19 , where c is the actual speed 
of sound, c_ is the isentropic speed of sound, a is the pipe radius 
in inches, and 5 is the kinematic viscosity in in2/sec. 
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to below 1000 cps. It is seen that relatively high accelerations can 
be obtained which if they occured at low enough frequencies or if 
they corresponded to fluid resonant frequencies would cause pressure 
signals large enough to be of great concern. From this Figure the 
pipe damping ratio can be seen to be about .04. 

Now the same system was excited with white noise input. Figs. 
2.13 and 2.14 show the pressure response for reservoir pressures of 
100 and 200 psia respectively. Notice that the two Figures are 
exactly the same except the first is everywhere between five and 
seven db less than the E N: This is exactly what would be expected 
as the induced pressure is proportional to density which for a gas is 
proportional to pressure. Thus when the pressure is doubled there 
should be a rise of six db in induced pressure noise, Fig. 2.15. 
These two Figures also show the fluid resonant frequencies at 60, 
200, 355, 455, 600, and 740 cps. Notice that the resonances at 455 
and 740 cps are much smaller than the others. This showed up in the 
theory also. The resonance at 600 cps is higher than it would 
normally be since although white noise was put into the vibration 
generator, the generator had resonances around 3/0 and 530 cps which 
show up throughout the graphs. The peak at 850 cps is due to the pipe 
longitudinal resonance. 

The several runs attempted with water filled pipes produced 


pressure levels approximately four times as high as for the 200 psia 


An explanation on reading db graphs is given in Fig. 2.15. 
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air for a given excitation. The increased density and speed of 
sound would predict level increases by a factor of ten. However, 
increased damping lowers the resonant peaks so that the factor of 
four seems reasonable with the theory. 

There is a misleading upward trend to all the graphs made by the 
analyser and recorder. This is because the filter bandwidth is 
always a constant percentage of the mid frequency being swept, in all 
these cases 6%. Thus the filter bandwidth doubles each octave, and in 


the case of white noise, twice as much signal is being averaged. 


oct 


2 
Ë 


This 3 db rise can be seen in Fig. 2.16 which is the random white 





Thus there is a 3 db rise per octave since i0 log = 10 log 2 = 3. 
noise input to the vibration generator. 
Fig. 2.1/7 shows the acceleration at the nonforced end of the 
pipe. The peak near 900 cps is due to the pipe longitudinal resonance, 
and the two peaks at 3/0 and 530 cps are the resonances of the shaker. 
In order to better understand these last several results, Fig. 2.14 
has been replotted in Fig. 2.18. In this Figure the 3 db rise has 
been eliminated and the pressure has been normalized by the input 
acceleration. Theoretical results are also plotted and agree 
favorably. This Figure shows the big drop in pressure after the 
first two resonances, and also shows that the resonances at 455 and 
740 cps are almost nil. The resonance at 600 cps is much higher than 
predicted as explained before. Because of the wideness of the filter 
at higher frequencies, the resonance at 8/75 cps is masked in the pipe 


resonance in the experimental results. 
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Before leaving longitudinal excitation, Fig. 2.19 is included to 
show that the fluid resonances are dependent on the total pipe length 
and not on the forcing fiow source positions. For this case, the 
distance between flow sources (pipe bends) was left the same as above 
while increasing the total distance between the reservoir and pressure 
transducer. It can be seen that the fluid resonances are at lower 
frequencies, 40, 156, 250, 350, 450, and 550 cps, corresponding to the 
greater length. As before the high peak at 550 cps is due to the 
Shaker resonance and the peak near 850 cps is once again the pipe 
resonance. 

Figs. 2.20 and 2.21 show the pressure response to random transverse 
excitation for clamped-clamped pipes of lengths 2/ and 23 inches 
respectively. The overall lengths from reservoir to pressure transducer 
are 46.1 and 42.1 inches, and in both cases there is zero mean flow. 
These pipes were not excited by inertial forces as the vibration 
generator was not large enough to move the whole system. Instead the 
pipes were forced directly at a point near one end. This method 
proved satisfactory. However, since the forcing was not symmetric as 
in the case of inertial loading, all the resonances are present and 
not just the odd numbered ones. 

Notice the difference between the two Figures. The first has many 
small resonances while the second has several large ones. This is 
because in the case of the 23 inch pipe, the fluid resonances and pipe 
bending resonances coincided. The two graphs are similar in that the 
first resonances are large while the higher ones fall off quickly. 


This is more noticable if the 3 db rise per octave is subtracted. This 
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fall off was predicted and is much more severe than for longitudinal 
excitation. 

Now looking at Fig. 2.20 in detail, the first peak at 44 cps is 
excited by the first bending resonance for a pipe of this length with 
clamped-pinned end conditions. Therefore, even though the pipe ends 
where clamped, torcing the pipe near one end excited this mode of 
bending also since the support moved torsionally. The peak at about 
86 cps is due to the double resonant frequency as predicted. The 
second and third resonant frequencies cause the peaks at 130 and 
260 cps respectively, but both their double frequencies are masked 
by other peaks. The one is masked by the chird resonant frequency 
and the other by the shaker resonance around 530 cps. The peak at 
5/ cps is excited by the tirst resonant frequency for clamped-clamped 
end conditions. its double frequency snows up near 120 cps. The 
fluid resonances are at 68, 220, 368, 515, 660, and 803 cps respectively. 
Only one, 517 cps, corresponds with any of the acceleration peaks. 

In Fig. 2.21 the fluid resonances can be observed at 80, 242, 
410, 585, 741, and 920 cps. The bending resonances for clamped- 
clamped end conditions are at 73, 216, 420, and 691 cps. Since the 
first several correspond aimost exactly with the fluid resonances, the 
peaks are seen to be relatively large. A small peak is seen at the 
first double resonant frequency also, about 152 cps. There is some 


motion of this beam in the clamped-pinned mode, but since these 


The pipe bending resonances were computed by directly using the 
formulas listed in Den Hartog [ 20 | : 
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frequencies do not correspond to any fluid resonant frequencies, the 
small peaks are masked out, 

Fig. 2.22 shovs the longitudinal accelerations at one support of 
this pipe. For the clamped-clamped mode the first resonance and its 
double frequency show at 73 and 158 cps respectively. Likewise for 
the clamped-pinned mode the first two resonances and their double 
frequencies are seen at 58, 111, 184, and 360 cps. The peak at 39 cps 
must be due to the first resonance of the pinned-pinned mode. 

In these last three Figures it is seen that the primary bending 
frequencies do excite pressures through longitudinal motions. Also, 
the double bending frequency accelerations are present, but they are 
of smailer magnitude and only show up for the first few resonances. 


Table 2.4 compares some of the experimental and theoretical results. 


TABLE 2.4 


Comparison of Experimental and Theoretical Ratios of Pipe End 


Longitudinal Acceleration to Platform Transverse Acceleration* 


(pipe is 1/4" copper tubing filled with air at 200 psia) 











Clamped-Clamped Clamped-Pinned 
g 8 8 8 
Cps Exp „end Theo end Cps Exp end Theo — 
base Sbase base base 
120 32 "za 98 „40 "Pa 
156 26 „62 Tı 23 ZB 





= 
The experimental datahad to be converted to g's per base acceleration 
since the pipe was not inertially forced. This was done by using 

Eq. (2.5-2) and a similar one for clamped-pinned end conditions. 
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This Table shows the theoretical results to be a little high 
because the end support stiffness in reality places some restraint 
on the motion. However, the results are of the right order of 
magnitude and show that the techniques of Section 2.5 could be used 
to determine an upper bound for the end accelerations. 

On the other hand the magnitude of the longitudinal end accelera- 
tions due to the primary bending resonances can not be determined at 
the present time. Considerably more experiments with various end 
conditions and base supports will have to be conducted in order to 
be able to predict this effect. The most that can be said now is that 
the first resonance is the most important and seems to induce a 
pressure peak of approximately the same magnitude as that induced by 
comparable longitudinal excitation. 

The last two Figures included are to show the combined effects 
of both longitudinal and transverse accelerations. The pipe was 
excited near one end at an angle of 45 degrees to its axis, giving 
both types of forcing. Fig. 2.23 shows the pressure response to have 
the general upward trend, on the recorder graph, of longitudinal 
excitation, with fluid resonances near 80, 242, and 400 cps. However, 
the bending resonances are also present. For the clamped-pinned mode 
the first resonance and its double frequency are at 59 and 121 cps 
respectively, but those for the clamped-clamped mode are masked behind 
peaks. They can be seen, however, in Fig. 2.24 at 77 and 154 cps 
along with the second bending resonance at 216 cps. These two Figures 
show that both longitudinal and transverse vibrations can excite fluid 


pressures at the same time, making the system susceptible to noise 
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regardless of the orientation of the system vith respect to the 
vibration field. 
Finally, Table 2.5 is included with a summary of the experimentally 


determined damping ratios. 


TABLE 2.5 


Experimental Damping Ratios 


Fluid-no mean flow Fluid-mean flow Pipe-long Pipe-trans 
(200 psia) (200 psia, 14 ft/sec) (1/4" copper tubing, air at 200 psia) 
Gpe öz Ges ÖL Cs с cps 2 
45 „03 42 „04 850 „04 73 „04 
80 203 155 „04 900 „04% 78 ә 
22 „04 197 „04 
390 „04 350 „04 


It can be seen that the values are fairly consistant, having 
a slight increase with frequency. Also, the fluid damping with no 
flow is smaller than with flow since the orifice adds resistance when 


it is in the system. 


2.9 SUMMARY 
In Chapter 2 it has been shown that piping systems containing 
fluid can excite large pressures when their Supporting environment is 
subjected to vibrations. Energy transfer to the fluid is by means of 
This value was corrected from a recorder output value of .07. 


Since the filter bandwidth is approximately 50 cps at this point 
it gives a misleading broadening of the peak which must be corrected. 
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flow sources caused by pipe boundary nonuniformities such as bends 
and area changes moving relative to the fluid. When subject to a 
white noise environment the pressure spectrum is composed of numerous 
resonant peaks over a wide band of frequencies, the peaks being due 
to the various resonances in the system, 

The piping supports can be forced longitudinally, transversely, 
or in combination. In all cases it is the longitudinal motion of the 
end nonuniformities that forces the tluid since the pipe transverse 
bending vibrations also cause longitudinal motion of its ends. These 
longitudinal end motions are excited at the bending resonant frequencies 
and at twice the bending resonant frequencies. 

The induced pressures due to longitudinal environmental vibrations 
were predicted and the experimental data gave excellent agreement with 
the theory. The excitation coupling to the twice bending resonant 
trequency end motion was first observed and then a theory developed 
which allows one to predict pressure levels with a fair degree of 
accuracy. The coupling effect at the primary bending resonant 
frequencies was not anticipated, but the theory for longitudinal 
vibrations can be applied once the coupled longitudinal motion has 
been measured. The experimental results show this coupling to be such 
as to induce pressures of the same order of magnitude as those induced 
by comparable longitudinal excitation. 

Since both longitudinal and transverse accelerations of the 
piping supports induce fluid pressures, the system is susceptible to 
noise regardless of its orientation to the vibration field. Motion 


components in more than one of the directions simply cause an additive 
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effect to the longitudinal forcing motion. 

The level of pressure induced for a given excitation was seen 
to be proportional to the effective speed of sound in the fluid and 
to the density which is proportional to the pressure for compressible 
fluids. For the particular system investigated, air at 200 psia, 
pressures of approximately .6 psi per g of environmental acceleration 
were observed for the first several resonant peaks. These pressures 
would have been even greater for fluids of higher density as was 
indicated by the scanty data obtained for water filled pipes. 
Reference [13] contains acceleration spectra for aircraft, rockets, 
and ships which can be used to estimate the levels of induced pressures 
that systems in actual operation can experience. 

Extremely high pressures can be excited when a fiuid resonance 
coincides with that of the piping system. For a more complicated 
system it is impossible to distinguish between these two types of 
resonances, but there are still absolute system resonant peaks that 
are higher than the others. This phenomenon can cause ә 
induced pressures at the higher frequencies even though the fluid 
resonant peak levels normally fail off sharply after the first two 
or three. 

Experimental results show the fluid resonant peaks to be fairly 
sharp, and both the fluid and piping damping ratios appear to be .04 
for the cases tested. 

The good agreement between the theoretical and experimental results 


indicates that the assumption made neglecting the viscous boundary 
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layer transfer from pipe to air is permissible over a wide band of 
frequencies. 

Reducing the production and transmission of these induced 
pressures is not simple. Use of a typical acoustic filter would not 
be satisfactory since the filter connections would only furnish more 
nonuniformities with which to generate more noise. The use of a 
compliant section at the valve end of the supply transmission line 
might, however, prove to be useful. Also, low pressure systems are 
less susceptible to the noise due to their low density, and perhaps 
should be exclusively used where very high accelerations will be 
encountered. Securing the piping to stiff structure at short 
intervals would also reduce the transverse vibrations and make the 
system susceptible to vibrations in only one axis, longitudinal. 
Lastly, a thorough vibration test of the empty piping system should 
be made since high induced pressures can be expected at any frequency 


where strong structural resonances exist. 
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Figure 2.1 Basic Fluid System 
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Figure 2.2 Fluid Control Volume 
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Figure 2,3 Pipe Element 
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Figure 2 4 Sample Problem 
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Figure 2.6 Shape of Damped Resonant Peaks 
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Figure 2./ Test Setup for Longitudinal Pipe Vibrations 





Figure 2.8 Test Setup for Transverse Pipe Vibrations 
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Figure 2.9 Experimental Arrangement 
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Figure 2.15 Explanation and Scale £or Reading db 


db = 10 log Se where q. is the reference quantity. Spectrum, 
O 
however, are in units of quantity squared. Therefore, 


2 
db = 10 log E = 20 log oa 
do o 


Graphs that are presented in db can be easily read by using the 
scale below. The correct quantity for a given db is obtained by 
multiplying the reference quantity by the appropriate factor in the 
scale. For units between decades, the decade factor is first multi- 
plied by the appropriate unit factor and then times the reference 


quantity. i.e. 23 db = 10x1.41 de 14.1 ët 


db qox db q x 

0 1.00 20 10.0 

1 5 30 31.6 

2 6 40 100 

3 1.41 50 316 

4 1.58 60 1000 

3 1.78 70 3160 

6 2.00 80 10000 
7 2.25 90 31600 
8 21 100 100000 
2 2,82 110 316000 


10 3216 120 1000000 
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3. CONTROL VALVE ANALYSIS 


3.1 INTRODUCTION 


In Chapter 2, fluid noise introduction through the environmental 
vibrations of a passive transmission line element was developed. In 
this Chapter an active element, a flapper-nozzle control valve, will 
be investigated while in the same vibrating environment. There are 
two transmission lines associated with the control valve, the supply 
line leading to the upstream orifice and the load line that transmits 
the signal to the load orifice. Noise can be introduced into the 
system in either or both of these. Induced relative motion of the valve 
Flapper also creates a noise problem and will be investigated along 
with the two other sources. The small dimensions of a control valve 
prohibit the generation of structural induced internal pressures since 
the resonant frequencies are so high as to make the pressures 
insignificant. 

A quasi-linear approach will be used whenever nonlinearities are 
encountered by linearizing the relation around the operating point. 
This method gives satisfactory results for small excursions of the 
variables. 

In the following sections the valve chamber and flapper will 
be analysed first. Then two different models, distributed and lumped, 
will be constructed for the load line. Finally, the experimental 


results will be presented and discussed. 
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3.2 MATHEMATICAL MODELS 


3.2.1 VALVE CHAMBER 


In a flapper-nozzle control valve the variable orifice is the 
curtain area between the nozzle and flapper. This orifice is located 
downstream of a fixed orifice so that any change in its area alters 
the mass flow rate out of the valve and causes the pressure in the 
chamber to vary, see Fig. 3.1. For the valve analysed, both the fixed 
and variable orifices have choked flow through them and behave as pure 
resistances. 


Writing the continuity equation for the valve chamber of Fig. 3.1 


gives 
w +w, -vW sa p V = О (3.2-1) 
4, dt c c ? 
where w = ... and v = (?22Nn"Dix,pe are the choked mass 


m. n 
| T VT 
O O 
flow rates through the upstream supply orifice and the variable nozzle 
orifice respectively. In these relations the N's are the orifice 
discharge coefficients, T. 1s the absolute temperature, and A. and 
TD x, are the areas of the respective orifices. Ҹу is the mass flow 


rate into the transmission line, and the valve chamber volume, Vos is 


taken to include that of the bellows. 





—::: Psi 
Let EE = be the initial supply steady 
\/ I s 
O 


state mass flow rate at the design operating point, where R. is the 
resistance of the orifice, and the subscript i is used to denote 


initial steady state condition. Then using the equation of state, 
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Pc 


də p RT 
О О 





, and the differential operator, D, Eq. (3.2-1) can be 


nondimensionalized to give 


хр, p p P ¿PX 


SE 5 





N = 8 + TD + tb = () (3.252) 
riPsi Psi Psi si Psi ri 
and for the steady state, 
= 7 
Poi É Psi (3.2-3) 
pP. 
si d 
In these equations 
Peco hee инни s Pes b 
NA : w.R T : b w.RT 4 
5 5 i oo ioon 


where R. is the gas constant and Ар is the area of the bellows. Тһе 


lengths &, and ei are defined in Fig. 3.1. 


b 


In order to get rid of the nonlinearity in the first term, 


Eq. (3.2-2) can be put into a quasi-linear form. Thus, 





Pei Ах, Ар, Aw, Ар. 
€—(r D 4 N) — + (TD + N) -— R ——— - — = 0 
A e x. De S : DE: 
si ri si si si 
Or Ap poe Ax Aw 
— _— pn) —- R — 
Pe Psi Psi : Xi S Poi 
— = —— w ә. (3.2-4) 
Psi TD + N 


In this equation it is seen that the valve chamber gives a first 
order lag to the pressure response. For optimum operation, the time 
constant T should be small and the signal to be transmitted introduced 
through the flapper displacement, Ах. However, with the valve in a 


vibrating environment, noise can be introduced by the transmission lines 
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in the form ot Ар. and Aw» and also as an additional induced 
flapper displacement. 

Ар. is the variation in pressure at the upstream orifice induced 
by the supply line vibrations. Since the orifice is choked, its 
impedence as seen by the line is a constant resistance and independent 
of the vaive characteristics. Therefore, Ap. can be independently 
determined by the methods of Chapter 2 and applied as a forcing 
function in Eq. (3.2-4). 

On the other hand the flapper motion, Ах, is very dependent 
on the rest of the control valve variables. This is also true of Aw, 
For although it is a transmission line induced mass flow, it cannot 
be independently determined since the impedence as Seen by the line 
is the valve chamber, This complicated impedence is wholly dependent 


on the other variables. Therefore, to completely describe the system, 


relations for the flapper and load line will have to. be developed. 


3.2.2 VALVE FLAPPER 


Consider now the behavior of the flapper. The force of the air 
jet acting on it can be found by applying the noninertial momentum 
* 
equation to the control volume in Fig. 3.2, and taking into account 


the unsteady terms which previously have been left out [23] . This gives 


The noninertial momentum equation is 


2 
z f v,(pv, dA) . 9 x dv P һ// v,(pdV) 
c.s. 50 E E Bo a с̧м So 


where the r stands for relative to the accelerating control volume. 
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(3.2-5) 
where fc is the force of the flapper on the fluid (or the fluid on 


the flapper), m. is the mass of air in the control volume, and x, is 
the motion of the platform. 

The above equation was obtained by noticing that the pressure 
forces cancel each other everywhere except over the nozzle area, A» 
Also in the first term on the right hand side, it is assumed that the 
air flowing outward through the curtain area has a mean axial relative 
velocity of one half the flapper relative velocity. In the last term, 
the integrand vanishes everywhere except in the volume between the 
flapper and nozzle. In order to evaluate this term a simple approxi- 
mate, one-dimensional flow solution can be integrated over this volume. 


This has been done in Appendix B with the results that 





Zë We = SSP AZ ) (3,926) 
dt 8, ot 3 2. 

CV. 

Also from the equation of state, 
Z 
p_v 
2 nn _ 2 

ün RT, o E —. 


vhere k is the ratio of specific heats, and M is the Mach number. 
Since the chamber area is much larger than that of the nozzle, the 


relation for stagnation pressure can be used to give 
k 


k-1 
x (1 4 k 2 1 M Si (3.2-8) 
c 


Le 


p 2 n > 





a 


Thus substituting Eqs. (3.2-6), (3.2-7), and (3.2-8) into (3.2-5), 


one obtains 





1 + kM x m ra ðX W 
f.=pA = - p + a O "JE. IAMS. 
f en k an g at? 6 ot g 

(1 + Lu A 
2 n 4 
2 —” Er | 
says A : (3.2-9) 


The term in brackets is merely a function of k and M. and thus a 


* * 
function of A , where An is the curtain area for choked flow. 


A 
n 
However, 
| * 
AW 2. 
A D 
n n 


Therefore this term can be replaced with an emperical expression that 


4x 
is a function o. . This was done by fitting a second order 
n 
relation to it, Fig. 3.3. The result, 
ix, 2 
70000 — 0” + 1 
D 
n 
ах 
gives excellent agreement up to кә „7. Virtually all flapper- 
n 


nozzle valve operation is well below this point. 
Now writing Eq. (3.2-9) in quasi-linear form after substituting 


in the above relation, the differential operator, and W gives 


Ax m 





Af. = A [.243x + 1| EAP EE ay 
f n c mero. g O 
ri O 
wN p. w.Nx_.DAp 
_ 4 JMC Dax ú 2 iri 06 | 32210) 
So Psi 80 Psi 
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the expression for the force the fluid exerts on the flapper, where 


The flapper is positioned by a torque motor and bellows feedback, 
Fig. 3.1. Because the torque motor mass, damping, and stiffness are 
considerably greater than those of the flapper, a lumped model for the 
flapper can be used giving the system a single degree of freedom. 
From the Figure it can be seen that an inertial force will act on the 
flapper if the platform experiences accelerations in the direction of 
the control valve principle axis. This will also be true for the torque 
motor rotor if it is not symmetrically balanced. 


The differential equation for the balanced condition is 


m х, пр 
“ə E sj 
n O 
b k 
1 Je 02 Pe Ze 
= z. (m+ 2 )D 25 (bęt ie z )x (3.2-11) 


where A is the signal torque applied to the torque motor and ky and 


k. are the stiffnesses of the bellows and spring respectively. Mes 


јӱз bę» b kę» and k. are the mass, inertia, damping, and stiffness 


t ° 


of the flapper and torque motor, and the lengths £, 4 Rs and Re 


b? 
are defined in Fig. 3.1. When this equation is written in the quasi- 


linear form, put into Eq. (3.2-10) and the whole expression nondimen- 


Sionalized by kx a> one obtains, 


p” 11 55: pa; | Ах, 
AF - AI = | WAG мк тмә ) —— 4*(1 k- Sox ——) — 
x Psi = Poi ri 
A Ap 
+ = oN = a(.243X b: 1) a> 2 (3.2-12) 
3 u A p 
n 5 41 
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In this Equation, 














2 
.. АР, ә mD x. š m AG o kn kt. 
2kx : g kx, kx i kt : 
EL ri 1 n 
= nP si _ ki z De SE P BOK. 
kx > [e m m ; 
i gmk e e e 
e e 2 5 
° o 
Jt DE = 
where m, = m, T m, - m ; be = bę + 2 , and k, = ke BR 2 


Eq. (3.2-12) shows that the flapper motion is dependent on the 
chamber pressure. It also shows that the pressure is a function of 
the inertia force and susceptible to noise introduction by it. 

Still one more equation is needed in addition to (3.2-3), (3.2-4), 
and (3.2-12) to completely describe the system. This is the relation 
for the load transmission line. Models for this line will be constructed 


in the next two sections. 


3.2.3 TRANSMISSION LINE, DISTRIBUTED PARAMETERS 

The transmission line leads from the valve chamber to the choked 
load orifice. Depending on the vibration field, it can be excited 
longitudinally, transversely, or in a combination. Regardless of the 
excitation or the piping supports, the methods of Chapter 2 can be 
used to obtain a relation for Aw, in terms of Ар. Here the pipe will 
be assumed rigid and forced, along with the valve, in a longitudinal 
direction with an acceleration of р”лу . 


Using the relative transmission matrix from Section 2.4 for the 


fluid, one obtains 
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- ( i 7 
Ар, cosh I Ap p) + Z sinh b Aw (L) oP ADsinh r Ay. 
and 65 
il | 
Aw, > sinh T Ap(L) + cosh T Aw (L) + р АО(совһГ – 1) Ay,» 


C 


where Aw, is the relative mass flow rate out of the chamber, Eq. (3.2-1). 


Since the load orifice is choked, w, (L) = pu, where R is the 


orifice resistance. Therefore, 


Z 
1/Z sinh? + 1/R cosh? (1-coshr--*sinhT)AG 
Aw, = ( met = ыы 
cosh! + p Sinhr D(cosh[ + R SinhT) 
and in the steady state, 
Mer e 
Qi R | 
Putting this last expression into Eq. (3.2-3) gives, 
ci m. : (3.2-15) 
si S 
N + R 


Eqs. (3.2-14) and (3.2-15) along with (3.2-4) and (3.2-12) 
completely define the system. Eq. (3.2-14) is for one specific type 
of load line and excitation while the other three Equations are general. 
The system behavior will be analysed using these Equations in Section 3.3. 

Before going on to a lumped model development for the transmission 
line, it is of interest to consider the impedence the line sees looking 
into the valve. Once this is known the load transmission line is 
completely determined and can be independently analysed. The impedence 
only has meaning when the valve is passive or AF, AI, and Ар. are zero. 

Ар, Ар, 


= — , combining Eqs. (3.2-4), (3.2-12), 


For this case, since Z € — 
V Aw Aw) 


and (3.2-15) gives 





- 74 - 








2 
D D 
7 ә +a + + b) 
n 
2. = - - AA F —Ə N (3.2-16) 
ni D TP + N D 
(rD+N) C—,+ ар b)+( R ) (c və d) 
“n N + — a 
R 
11 N 1 2 
= ә — = X = a = 
where a = 27 + en R. , D 1 + k „486a R .. 3 nN, 
N+ — N+ e" 
Zn 
and d = a(.243X + + 1). 
5 


From Eq. (3.2-16), the impedence and real and imaginary parts are 
plotted in Fig. 3.4 for the particular control valve used in the 
experimental work. Also included is the impedence for when there is 
no bellows feedback, or Einer > A = 0. It is seen that feedback 
emphasizes the flapper resonance, but that it decreases the impedence 
at lower frequencies. 

It is to be noted that the real part of the impedence is positive 
up to the resonant frequency where it becomes negative, and the 
imaginary part is negative for all frequencies. Using this figure 
along with the known orifice resistance, Eqs. (3.2-13) can be used to 


solve the load line pressure and mass flow rate at any frequency. 
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MODELLING TECHNIQUES 
In addition to the distributed model for the line, it is also 


desirous to construct a lumped model which allows the system to be 
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analysed on the analog computer. The advantage of using the computer 
is that each of the parameters can be easily changed and its effect 
on the rest of the system readily seen. 

The line is modelled as lumped inertias, resistances, and 
compliances. The inertias and resistances are combined into series 
impedence elements and the compliances constitute shunt admittance 
elements, see Fig. 3.5. 

Causality must be taken into account when constructing the model 
so that the problem can be solved on the analog computer without using 
differentiation which causes difficulties. The specified variable, 
pressure or mass flow rate, at each end of the line determines what 
the nearest lumped element must be. Pressure requires a series 
impedence while mass flow rate requires a shunt admittance. Also, the 
number of line resonances a person desires to model determines the 
number of inertia elements that must be used. Thus since the trans- 
mission line in this investigation starts with a known pressure, Ар. 
and terminates with an induced flow source due to the motion of the 
line, and since it is desired to model the first two resonances of the 
line, the model in Fig. 3.5 is appropriate. 


The transmission matrix for the first half of the line is 


ae 


АР, 1 5 1 0 AP, 
= TE 
Aw» 0 ilı 2 1 De | 


and that for the second half is 
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бр. 1 = 1 O| |Ap(L) 
Aw 0 1 e 1 Aw(L) 


In relative coordinates, since Aw = Aw + P ADAY_» these become 


YLZL ZL ZL 
Ss e 4 2 Pr 2 Һу, 
= + 
Bore 2L 1 OF mn 0 
2 
and 
YLZE AL ZE 
Ар. 1 + ; > Ap(L) > p ADAy_ 
= + 
Aw = 1 SEND 0 
rm 2 R 
: Ap(L) . 
respectively, when the substitution Aw (L) Er is made. 


Now extending the results in Section 2.4, the series impedence 
becomes 


ZL = ~~ =R 


| 
N 

| 
+ 
A 


Ge (3.2-17) 


where T = -. the transmission time for a vave to travel the 
length of the pipe, and Ry is the total frictional resistance. Also, 


the shunt admittance becomes 
8 ADL 


2 
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(3.2-18) 
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Putting Eqs. (3.2-17) and (3.2-18) into the transmission matrices 


then gives 2 9 a 2 





T D RTD Z T D R Ar p AD Ay 
Se e 2 — m5 Ap ge O 
c 4 62. 2 2 m 2 
E r + 
Aw — 1 Aw 0 
(3.2-19) 
and 
* 
Zei R Tp D Ry АД 
Ар. mr “Wi 575:7—UCUOC(7:— —5 
| E | 
= |T D | | 
5... 1 5. 0 | 
m | 22 R = 
C a 


For the steady state, Eqs. (3.2-19) give 
x Poi 


7:1: R, 


which when combined with Fq. (3.2-3) becomes, 





1600001 ары" (3,2-20) 
51775 и 
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Although Eqs. (3.2-19) could be cascaded to give one overall 
transmission matrix, this should not be done for the same reason as 
mentioned previously, namely the resulting equations would not be 
able to be solved on the analog computer without the use of differenta- 
tion. Therefore, regardless of the number of lumped sections the line 


is divided into, each section should be handled separately. 
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* | 
The resistance is not included in the forcing term since only 
the flow relative to the pipe causes a loss of pressure. 
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In order to determine Ro for use in the above equations, it is 


necessary to use a quasi-linear approach since the frictional resis- 


tance is nonlinear. The pressure loss equation is 


2 
4fLp V 


m 28D, 
where f is the friction factor, Po is the mean fluid density, and 
D, is the pipe diameter. Now putting in mass flow rate gives 
Е. AfLo w“ 
28 D,p A’ 
where A is the flow area. 


Putting this into the quasi-linear form it becomes 


8fLp w, Aw 
77 22 
28 .О,о, А 
Therefore, 
A 4fLw, 
bern 777 (3.2-21) 
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Using Eqs. (3.2-4), (3.2-12), (3.2-19), and (3.2-20) it is now 
possible to set up the problem on the analog computer. An operational 
block diagram for this setup is contained in Fig. 3.6. Since every- 
thing is nondimensionalized, including time, scaling is not too diffi- 
cult and is indicated on the Figure. The various inputs are also 
marked as well as the system variables that are of interest. It is 
noted that the load line induced noise input, Ay > enters in at two 
different points due to the load line being modelled in two parts. 

This analog model now allows one to vary the input and system 


parameters, and to investigate their effect on the system variables. 
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Results thus obtained are discussed along vith the other analytic 


and experimental results in the folloving section. 


3.3 EXPERIMENTAL WORK AND DISCUSSION OF RESULTS 

A small flapper-nozzle control valve was built in order to 
verify experimentally the analytical results. It was equipped with 
a removable bellows feedback system, Fig. 3.7. The flapper displace- 
ment could be measured by means of a capacitance pickup, and the 
chamber pressure and supply line pressure just upstream of the ori- 
fice were obtained with a piezoelectric transducer. The flapper 
displacement was controlled by means of a torque motor supplied 
from a servo amplifier. A view of the general arrangement and instru- 
mentation are shown in Fig. 3.8, and further discussion is contained 
in Appendix C. 

Since a vibration generator large enough to shake the whole 
system was not available, it was necessary to simulate the noise 
inputs. The pressures in the supply line were generated by supplying 
the valve through the load line of another flapper-nozzle valve 
operating at a higher pressure. In this manner pressures at any 
frequency could be imposed on top of the mean supply pressure of 
100 psia, simulating resonant noise peaks. 

The behavior of the control valve transfer function of, Ар. 
when subjected to supply line noise, AP» is shown in the frequency 
response curves of Fig. 3.9. The Figure contains an experimentally 
obtained curve for the valve operating with bellows feedback as well 


as theoretical curves with and without feedback. 
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The theoretical results were obtained by combining Eqs. (3.2-4), 
(3.2-12), (3.2-14), and (3.2-15). After putting in the appropriate 
values for the valve used, this distributed model equation becomes 


adp. + .32bAG + +Ol2acAG, 


MPe “ae + 6.75b + 1.82ad | re 


where Ap. is in psi, Ар. is the supply line pressure noise in psi, 
AG, is the acceleration in g's experiences by the valve in the direction 
of its principle axis, and AG, is the acceleration of the load line 


in g's. Any one of the three noise transfer functions can be ob- 


tained by setting the other two noise inputs to zero. In this 





equation, 
W | 
a = 4000 1.12: "=" 
u) 2 
n W 
n 
b = .02 — i + .22, 
W 
n 
“a W W | 
c = — | (cos 1100 7 1)i + .23d sin 1100. >» 
O S + 4,40 i si — 
55 1100 E 100 
dE W W i 
cos 1100 * 5257 1100 
e = 7.36 — í + .50, 
n 
and 


мә 2520 radians per second. 


When there is no bellow feedback, the denominator of Eq. (3.3-1) 


becomes ae + 1.30b + 1.82ad, 


2 

.35 — i + 1.00 - —, 
w w 

n n 


and a 


while b = .22, 
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In the Figure, the transfer function is seen to behave as a 
first order lag cascaded with a second order system. This is due to 
the valve chamber and flapper respectively. Feedback increases the 
effect of the flapper resonance, but it decreases the pressure at 
lower frequencies. The experimental curve is not continued past 
700 cps because the response of the first (noise producing) flapper - 
nozZle valve drops off to zero at this frequency, and no noise could 
be introduced to the test valve. 

The system appears to have slightly more damping than was 
estimated, and therefore the experimental peak is not as high as 
predicted, and it is at a slightly lower frequency. From the Figure, 
the damping is calculated to be approximately .08 and .10 for the 
theoretical and experimental results respectively. It is to be noted 
that the flow noise in the valve created a region of uncertainty in 
the experimental results having a width of approximately .02 on the 
Figure. 

It appears that if the flapper-torque motor arrangement had less 
damping and if the bellows feedback were increased, then the noise 
transfer function would become much larger near the flapper resonance. 
In order to study the effect of varying these two parameters, analog 
computer solutions were obtained using the lumped model formulated in 
Section 3.2.4, 

It was found that the noise signal could be increased as expected. 
However of much more interest, it was also found that the total noise 


power could be amplified in the valve at frequencies near the flapper 
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resonance. That is to say that the noise pover associated vith the 
bellovs, force times velocity, vas greater than the noise pover 
entering the valve, pressure times flow. This amplification can take 
place because some of the steady power flowing into the valve is 
converted into noise power. This is a danger of noise in an active 
element and differs from what can happen in a passive element. Ina 
passive element high noise signals can be obtained through resonances 
as in the transmissian lines studied, but the total noise power cannot 
increase. It is of interest that the increased power was measured in 
the bellows which is a direct analogy to a load piston operating with 
the control valve. Thus it is seen that the noise car be amplified 
into actually displacing a load. Also, increased motion of the 
flapper due to the feedback of the supply line noise can cause the 
flapper to strike the nozzle. This flapper motion saturation makes 
the valve inoperative to command signals. 

In terms of numbers, the threshold of amplification was observed 
when damping was decreased by a factor of .7 and feedback was increased 
by a factor of 1.7. Another combination was found when damping was 
decreased by a half and feedback was increased by a factor of 1.2. 
These factors decreased the flapper system damping ratios to .25 and 
„17 respectively for the two cases. These are not excessively low 
values and might well be encountered in many systems. These results 
could not be checked experimentally using the test valve because 
the feedback and torque motor damping had fixed values. When the torque 
motor was removed to decrease the amount of damping, the system became 


unstable. Since the feedback could not be increased, it also was not 
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possible to make the flapper motion saturate. The maximum flapper 
motion that could be excited vas about .0015 inches vhile the steady 
state nozzle clearance vas .005 inches. 

In order to simulate inertial excitation to the flapper, a vhite 
noise signal vas added into the servo amplifier input. This shoved 
the first resonant peak as expected, but it also shoved smaller 
peaks at the equalizing spring resonances and at the second resonance 
of the flapper. Fig. 3.10 is the frequency response transfer function 
of the chamber pressure per g of induced flapper excitation, The 
characteristics are very much like those discussed in coniunction 
vith Fig. 3.9. There is, hovever, a small peak at 300 cps in the 
experimental result due to the equalizing spring, and the pressures at 
the lower frequencies are larger than predicted. Increasing feedback 
and decreasing damping would have the same effects with regards to 
amplifying this induced noise as was discussed with supply line noise. 
It is also possible to get flapper motion saturation. 

Since there was no driver available that would work against the 
pressure in the load line, it was not possible to simulate load line 
induced vibrations. The line was hooked up to the vibration generator, 
but this did not induce signals large enough to be recorded. Fig. 3.11 
shows that the theory predicts that these pressures will be very small. 
Even for 100 g's of excitation the induced pressure would be less than 
one psi even at the flapper resonant frequency. The smallness of the 
induced pressure is due to the capacitive effect of the valve chamber 
volume, and to the low pressure of the load line, approximately 35 psia, 


causing the air density to be quite low. A valve operating at higher 
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pressures would be much more susceptible to this type of noise. 

The above results have shown that the control valve is susceptible 
to noise from the supply line and inertial forcing of the flapper. 
The noise amplitude is greatest at low frequencies and near the flapper 
resonant frequency. When the flapper system damping is low enough, 
bellows force feedback can cause amplification of the noise power 
to the valve and saturation of the flapper motion. 

External vibrations of the load line do cause resonant fre- 
quency peaks, but because of the reduced pressure of this line, it is 


only of concern invalves with very high supply pressures. 


3.4 PROCEDURE FOR TESTING CONTROL VALVE SUSCEPTIBILITY TO NOISE 
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Of considerable interest to the person who must specify a parti- 
cular control valve for a given job, is whether or not the valve will 
be highly susceptible to vibration induced noise. Therefore, it should 
become standard procedure to test any such control valve before it is 
assigned for use. Since there are three ways noise can be induced 
into the system, there are primarily three tests. 

The first test deals with noise introduced through the supply 
line. For this test, the supply line noise can be generated in the 
same manner used in this theses. Supply the valve through the load 
line of a second control valve operating at a higher pressure. In this 
way sinusoidal noise signals can be imposed on the mean supply pressure 
for any frequency. The frequency response for LG (w) | should then 


be obtained as well as the response for 1— il, making sure to test 
S 
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at the maximum feedback that will be used in practice. If either of 
these two responses is excessively large, the control valve should not 
be used where environmental vibrations could be a problem. 

The second test is for vibration induced flapper motion. If a 
large vibration shaker is available the valve can be vibrated directly 
and EE [anal | recorded. Otherwise the frequency response 
for both of these can be obtained by applying a signal to the torque 


AF 


motor as was done in this thesis. In this case G= — where F 
f 


is the torque applied to the motor, m, is the equivalent forced mass 


f 
of the flapper, and 4 is the lever arm of the rotor. If either of the 
responses is too large, an effort can be made to balance the rotor and 
flapper assembly. If one is fortunate enough to have the large shaker, 
any reduction in response due to balancing will show up directly. 
However, if one is simulating the vibration, balancing will reduce the 
value of M used in calculating AG. If M is reduced enough, AG will 
be excessive in order to force a response, and the control valve will 
be safe for use. If balancing is not successful the valve should not 
be used. 

The third test is for load line induced pressures. If the control 
valve chamber does not operate above several hundred psia, this type 
of noise should not be a problem. A test can be made with a small 
vibration generator used to vibrate the end of a flexible line in a 
longitudinal direction and the response, Fra) | recorded. The 
amplitude of this response will dictate whether the valve is safe or 


not. 


One additional test can be conducted with the fluid transmission 
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lines. A thorough vibration test of the empty piping system should 

be made since high induced pressures can be expected at any frequency 
where strong structural resonances exist. If any of these resonances 
are near the flapper resonant frequency of the proposed control valve, 


a different valve should be used. 
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Figure 3.1 Control Valve Arrangement 
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Figure 3.2 Flapper-Nozzle Control Volume 
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Figure 3.3 Nozzle Relation Approximation 
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Figure 3.5 Lumped Parameter Model of Load Transmission Line 
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Figure 3.6 Analog Computer Operational Block Diagram 
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Figure 3.7 Flapper-Nozzle Control Valve 
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Figure 3.8 Experimental Setup 
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4. CONCLUSIONS 


4.3 VIBRATION INDUCED TRANSMISSION LINE NOISE 
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From the results of Chapter 2 it is concluded that piping 
subject to white noise environmental vibrations can excite resonant 
pressure peaks in the contained fluid over a wide band of frequencies. 
The pressures are caused by induced fluid flow at pipe boundary 
nonuniformities (bends and area changes) due to the relative motion 
of the piping. The peaks occur at the various fluid and piping 
resonant frequencies. 

Either longitudinal or transverse piping support vibrations, 
or a combination, will cause induced pressures. In all cases it is 
the longitudinal motion at the end nonuniformities that forces the 
fluid. This is true even for transverse vibrations since the trans- 
verse pipe motion exciteS longitudinal vibrations of the ends at the 
pipe bending resonant frequencies and at double the resonant 
frequencies. Since pressures are induced for either type of vibra- 
tion, a fluid piping system is susceptible to noise generation 
regardless of its orientation to the vibration field. 

The level of induced pressure is proportional to the fluid 
density, the effective speed of sound in the fluid, and the amount 
of excitation. Resonant peak levels of .6 psi per g of longi- 
tudinal acceleration can be obtained at the lower frequencies for 
air at 200 psia mean pressure in one quarter inch copper tubing. 


Preliminary results show that correspondingly higher pressures can 
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be obtained with liquids due to their higher density. The 
pressures generated by the first several pipe resonant frequencies 
when subjected to transverse vibrations are comparable to the 
longitudinally excited case. However, levels due to transverse 
double resonant frequency excitation are less by about a factor 

of two. 

Extremely high pressures can be excited when a fluid resonance 
coincides with that of the piping system. For a more complicated 
system it is impossible to distinguish between these two types of 
resonances, but there are still corresponding absolute system 


resonances that are higher than the others. 


4.2 CONTROL VALVE 


For the one type of active Zlece studied while in a 
vibrating environment, noise introduction into the valve takes 
place through the supply transmission line, induced motion of the 
flapper assembly, and for high pressure systems, through vibrations 
of the transmission line to the load. Since the valve transfer 
characteristic is that of a first order lag cascaded with a second 
order system, due to the capacitive chamber and flapper resonance 
respectively, the noise transfer functions are largest at the low 
frequencies up to the break frequency and again near the flapper 


resonant frequency. 


A Pneumatic nozzle-flapper control valve. 
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Force feedback in the way of a bellows arrangement reduces the 
pressure levels at the lower frequencies. However, it also intensi- 
fies those near the flapper resonance. By increasing the feedback 
enough, it is possible to get amplification of the noise power for 
frequencies near the flapper resonance. That is to say that in any 
active element it is possible to convert some of the steady flow 
power entering the valve into noise power. This phenomenon does 
not occur in a passive element where signal amplification can take 


place, but not power amplification. 
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5. RECOMMENDATIONS 
5.1 FUTURE WORK 


Since the results of this study have shown that significant 
vibration induced pressures can be generated in a fluid control 
system, further work in this field is warranted. The most fruitful 
areas for future work are the following. First, it is recommended 
that extensive experimental work with transversely vibrating pipes 
be conducted in order to better understand the structural coupling 
between longitudinal and transverse motion. Data should be taken 
for various pipe end constraints and different support geometries. 
Secondly, different types of active elements should be analysed 
while operating in a vibrating environment. This should include pure 
fluid elements such as bistable and jet amplifiers, as well as systems 
with several stages of feedback and amplification. The results should 
be compared along with those for the flapper-nozzle valve to determine 
which elements are least susceptible to vibration induced noise. 
Further experimental work should also be conducted with the flapper- 
nozzle valve to verify that noise power amplification takes place. 
Lastly, it is recommended that parallel work be conducted for a 
hydraulic system. Preliminary work has shown that a liquid system 


will be even more susceptible to vibration induced noises. 


5.2 APPLICATION 
Information learned in this thesis can be applied in order to 
reduce a systems susceptibility to noise. Thus, it is recommended that 


all fluid equipment that will be exposed to strong environmental 
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vibrations be thoroughly tested before use. The system should be 
shock mounted if possible. Otherwise it should be mounted on stiff 
Supports so that all resonant frequencies will be high. The 
transmission lines should also be secured at short intervals in order 
to decrease transverse vibrations. The installed system, empty 

of fluid, should then be vibration tested. High induced pressures 
can be expected at any frequency where strong structural resonances 
exist. Lastly, the fluid amplifiers should be carefully tested in 
order to pick those with good transmission characteristics that will 


not cause noise amplification. 
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APPENDIX 
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APPENDIX A 
TRANSMISSION MATRICES 
A.1 TRANSMISSION MATRIX DEVELOPMENT 


* 
A model for a uniform, one-dimensional, distributed system is, 
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where Z is the series impedance per unit length, and Y is the shunt 
ə 
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Using pressure, p, and mass flov rate, v, as variables (force and 


admittance. Either or both can contain the operation D € 


displacement for pipes) the relations for a micro-element become, 


dp _ _ 
AS Zw 
E a (A.1-1) 
"w 
or 
p(x) > 0 z p(x) 
d_ — 
dx u 
w(x) Y 0 w(x) 


* 
References [15 | and [21] give extensive developments containing the 


material covered in section A.1. 
kk 
When working with mechanical and fluid systems, this form of the 


equations is obtained when the pressure or force variable is taken as 
compressive and the mass flow rate or displacement variable is taken 


in the direction of x. i.e. Ss 
=> 
Fes Ÿ TT 7 
PA wi 







ə — = 
ә 
ke e, 
– 
| 
' 
a = 
4 ma ün ә . 55 
le — © 211... au q uv © — 
O OOO 85 — (| > ue Ы 


T mms en © ua» 
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Ordinary differential equations are written instead of partial 
because time is contained in the operator D in the system parameters. 


Writing the above matrix equation in compressed form and using 


matrix algebra, one finds that 


ad (x) = Фи 


dx 


$ (xə) 


Jx) 


H H 
> e 
о о 


and 


l, = ~ D (x). 


Therefore 


Ас x) 
2 1 
S (x) =e D (x,) 655 
vhere Co is the integration constant matrix. 
* 
The transmission matrix can be expanded into, 
_AG- хј) _ cosh y(x,- x) Z sinh y(x,- хј) 


1 
2 sinh y(x,- хј) cosh y(x,- Xy) 


where Z = İZ/y is defined as the characteristic impedence and 


у = \ YZ as the propagation operator. 


The final expression becomes, 


p(x,) cosh y(x,- x,) 2. sinh y(x,- x.) p(x,) 
5 Z= sinh y(x,- x4) cosh y(x,- хј) и(х,) 
с : 
(A.1-3) 


This can be done by using either the Cayley Hamilton or Sylvester 
theorems [22] : 
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This transmission matrix is applicable to any uniform, one-dimensional, 
linear, distributed system of two variables. When working with piping, 


p and w are replaced by f, force, and 6, displacement, respectively. 


A.2 TRANSMISSION MATRIX USING RELATIVE FLOW OR DISPLACEMENT VARIABLES 


——npı.cıə ë “rr rt atus qm“ ee ee eee 


There are times when working with a system mounted on a vibrating 
platform when it is more convenient to use flow or displacement variables 
: + : : 
relative to the platform motion. Thus for distributed fluid and 


mechanical systems the relations have the form 


DEE SZW 

dx r o 

dw, 

Кур (А.2-1) 


Since w(x) = w, (x) + Wo: The r and o stand for relative to the 
platform and platform respectively. The platform is rigid so that 
its motion is independent of position. 


Using matrix algebra one obtains from 


p (x) 0 Z p (x) ~ Zw 


d_ 
dx 
w, (x) Y 0 w (x) 0 


that 
> =- Å 5 Gi +C 
-A x cl 
D (x,) “CAC 


D Gp) KI C.+AC 


" ar = 
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and 
A\x 
(se ` Dr s" À > A | 


b 6) - dasi” (x5) + 1 кә Ы AC 


Therefore 


(A.2-2) 
Now since 0 1/Y - Zw 
-1 O 
A - and C = 
1/Z 0 0 
the final expression becomes, 
p(x,) cosh Y (XX Z sinh Y (XX р(х,) 
= + 
w,(x,) 1/Z sinh YÜx,-x)) cosh Y(x,-x,) v (xə) 
2. sinh Y (x57x)w_ 


[cosh Y(x2-x,)-1] 5 


It is to be noted that Eq. (A.2-3) can be obtained more directly by 


ә 


substituting the expression w(x) = w(x) + Wo into Eq. (A.1-3). This 
method is not recommended, however, for more complicated systems 
involving matrices of higher order. In those cases, application of 


Eq. (A.2-2) will be more profitable. 


A.3 TRANSMISSION MATRIX WITH TRANSPORT 


* 
By introducing wave-scattering variables one can find a good 


approximate form of the transmission matrix when there is transport 


*For a comprehensive presentation of wave-scattering variables see 
Paynter [15] , pp. 268-271, and Brown [21] : 
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present. This would be useful for cases of mean flow in fluid lines. 
Wave-scattering variables can be obtained from the previous expression 


45, = -AS by obtaining the canonical form of the differential 
-1 


equation by diagonalizing A\ . Thus, premultiplying by A gives 


Sek bech 
OS et annt 5. 
Defining с as E = "TİŞ 


Sl 
leaves | A TN as the desired diagonal matrix. 


Therefore, for a two variable system, 
wo E 
c c 
M == 
E LN = IN ә 
Le Ce 


) 


"Fez 


Mrs Më 


2 

and 

u(x,) N 457321 0 u(x,) 

z (A.3-1) 

V(x) 0 .... v(x.) 

since ; 
ü A (x -x / 

D (x) = M” e = MS &,) 

and 


u(x) ł -1 
vol 26 - PV Dw. 
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For a wavelike transmission line, i.e. fluid line, Eq. (A.3-1) 


becomes, 


u(x,) e 0 u(x») 


v(x) 0 e v(x.) (A .3-2) 


where T is the wave propagation time between points 1 and 2 and is 
Хд-Х 





equal to where c is some characteristic velocity of the wave 


propagation. When one looks closer it is seen that the scattering 
matrix contains only a pure time delay and a pure advance operator. 
Thus writing Eq. (A.3-2) in the time domain it becomes, 

u,(t) = u, (t-T) 


v) (t) = v (t-T), 


This means that u is a rightward traveling wave and v is a leftward 


traveling wave 


I: A 
> тг 
3 > 


as u, is u, delayed by T, and v 


1 1 55 v, delayed by T. 


Now if one introduces a mean transport velocity V, 


575: 
a > 
Ki V2 


the propagation times will be different for the two directions. For 
x,-x 





the rightvard traveling vave, T. = 1) ә 
XT Xi 


wave, T, = — °: Now Eq. (A.3-2) becomes 


and for the leftward traveling 
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u(x,) e 0 u(x,) 


v(x) 0 e v(x,) | 


and transforming back to the $ variables it becomes, 


ER VGA LAD ING V ` 
p (xj) 1 ze a e m 0 Ze ae p (x) 
22, 2 ИЛ Ar” "diir, AN 
w(x) Z - VZ 0 e Z c| |w(x,) 
E c c 2 


Multiplying gives 


p(x,) e +e Ze a ) r(x.) 


N |= 


w(x,) RZ - e ) e +e w(x) 


Now one can make the approximations that 





w ED ER x V 
242 (los 
CEF 
E 
and 
22 25 La TA 
2—:-.:mx v 
re H = uz 


since = << l for all cases where the above analysis holds. 


Substituting these relations back into the expression gives, 


p(x,) cosh y(x,-x, Z sinh y(x,-x] p(x,) 
-V(x_-x,)D 
25 
RÉ pe 
c 
w(x.) 1/Z sinh Y(x,—x)) cosh Y(x,- x) w(x,) 


(A.3-3) 
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Thus it can be seen that transport merely introduces a multi- 
plicative factor and that Eq. (A.3-3) approaches Eq. (A.1-3) as 
NB" 0. 

When the equations are written in terms of a relative flow or 


displacement variables one obtains, 


rx - = : = 
p (x) V( 2 x,)D cosh Y (x, Xy Z sinh Y (x, Xy 
c2 
= e z 
w, (x) 1/2 sinh Y (497%, cosh Y(x,-x) 


2 sinh TG xD 
+ V(x,-x,)D (A.3-4) 


| B c 
cosh y (x,=x} ) e w | 


In this case one of the internal terms is altered slightly. Yet this 


also approaches its corresponding no transport equation, Eq. (A.2-3), 


as V —0, 


р(х,) 


w, (x) 
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APPENDIX B 


UNSTEADY MOMENTUM TERM EVALUATION 


B.1 APPROXIMATE SOLUTION 


A simple one-dimensional flow solution can be integrated over 
the volume between the nozzle and flapper in order to obtain an 


approximate value for the integral, 


V 
d xr” 
[=a Í dV. 
öt g 
ev. o 


This is done by assuming that the one-dimensional flov lines are 
alvays perpendicular to the slant areas of the right circular and 
truncated cones into which the volume can be divided, Fig. B.l. 

The integration of this volume is done in two parts. The first 
includes the circular cone volume elements starting at the nozzle 
opening, and the second consists of the truncated cone elements 
which end at the curtain area. Dividing the problem this way gives 
an exact integration of the volune for the first part and a close 
approximation for the second. 

Beginning with the cones, let A = TSR be their slant areas 
where R is the radius and s is the slant distance. Then an element 


of volume is 
A sd0 
S 
3 





dV = 


where 5 is the factor for a cone, and 0 is defined in the Figure. 


It is also seen that 










—_ © emm 9 4 © <a 
> © -5£5—.::21-. 


* 
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R 
* cos 8 
w 
Thus, since p = ae , these expressions can be put into the first 
rs 


part of the integral to give 











8 
O 
ә Ур ә v,cosÚ ҹу ALR dé 
— dy = — —T OTT. 
ot g dt g v_A 3 cosé 
О) O O Y s 
O 
da 
w R dü 
с ERA 
= 9 
öt 3 Be 
O 
-1 "e 
where E = tan R ` Integrating then gives 
© -l 
c” wR dé _ SES w tan Ze 
öt 3 8, 3 ot Bo R 
O 
JD SEM 
3 ot Bo 
Er -l Xr Er 
since — is small and tan TR R : 


For the truncated cone elements, A. = ms(Rtr) is the slant area 
and r is the radius defined in the Figure. An element of volume is 


then once again 








A sdé 
dV = 3 
x, 
where s = SE 


Thus the remaining part of the integral becomes, 


LÉ 
2 


3 Ver? IV ә v cost Ss A х, 
öt g "är g v_ A_ sind ES 
O) O 9 O r s 


O 





SEIN 








au 
2 w x cot9 
7 de 
öt 2 i 
o о 
о 
Integrating, one obtains, 
m 
2 w x cote Wy r gn sind 
oe Eet, e EE 
öt 8 3 8, öt O 
o 
SE rr 
5. Bo \ x2r? 


x 
Since — is usually about .l to .3, this expression can be approxi- 
13 Hee Ху 
mated by 7 — ( 1546500707 1: 
3 öt 8, 





since the relative flow through any slant area is equal to the flow 


into the nozzle, or WWL 
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Figure B.1 Elemental Volumes 
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APPENDIX C 


COMMENTS ON EXPERIMENTAL PROCEDURES AND INSTRUMENTATION 


SSS ene me eee ee ee ee — –Д ҝЖУН— ` 


| s a i ey pe eee 


In Chapter 2, the physical set up for experimentally determining 
the induced pressures in vibrating pipes was described and figures 
provided. This Appendix contains some additional comments. 

While taking data it was necessary to secure the pressure trans- 
ducer to a nonvibrating support due to the sensitivity of the trans- 
ducer to vibrations. This was accomplished by having the copper tube 
terminate in a flexible section of poly-flow tubing which could then 
be secured. The other end of the piping was attached to the constant 
pressure reservoir in the same manner. 

The vibration generator was made from a 35 watt speaker unit 
with the cone removed. Extra care was necessary with alignment in 
order to keep the coil from rubbing while it was in operation. This 
speaker unit made a very satisfactory vibration shaker that responded 
well to input signals from both a sinusoidal oscillator and a random 
noise generator. These input signals were preamplified by a regular 
high fidelity amplifier before being fed to the shaker. It was 
necessary to monitor the current going to the shaker coil in order 
to keep from burning it out. 

The mean test pressure level was maintained by use of a regulator 
reducer that was fed from a high pressure air line. Any fluctuations 
were damped out by the constant pressure reservoir. A globe valve was 


installed upstream of the reservoir in order to bleed down the system. 


= ын к—лкжһ”/ — 2: = 
= 
anə E) = 
a” onı = 
@ :5 
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The pressures and accelerations were picked up with Kistler, 
Model 601-A and Endevco, Model 2215 piezoelectric transducers 
respectively. The signals then went through buffer amplifiers to the 
oscilloscope and analyser and recorder. The frequency analyser and 
level recorder were Bruel and Kjaer, Type 2107 and Type 2305 
respectively. This equipment automatically swept through the frequency 
range with a 6% filter bandwidth and recorded the mean square levels. 


Averaging time was controlled by the paper feed rate in the recorder. 


C.2 CONTROL VALVE NOISE EXPERIMENTATION 


The flapper-nozzle valve built for experimental purposes had a 
chamber volume of .1 cubic inches. The nozzle diameter was .116 
inches while those of the supply and load orifices were .086 and .116 
inches respectively. The valve was designed to operate at 100 psia 
mean supply pressure and with a steady state flapper displacement of 
-005 inches. 

The flapper displacement was controlled by a Midwestern, Type 
11 F torque motor that operated from a 300 volt d.c. American Measure- 
ment and Control servo amplifier, Model 652, The flapper displacement 
was measured by using a Decker Unit, Type 902-1. This capacitive 
pickup had the advantage of not imposing any unwanted force on the 
flapper. 

Pressure measurements were obtained with the same piezoelectric 
transducers previously mentioned. Readings could be taken in the valve 


chamber and just upstream of the supply orifice. The supply air was 





= 1 = 


provided through another flapper-nozzle valve whose load was the 
test valve. In order to maintain this mean supply pressure at 
100 psia, it was necessary to operate the first valve supply pressure 
at 467 psia due to the large mass flow rate required by the test 
valve. 

Data analysis was accomplished with the same equipment mentioned 


in the previous section. 
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